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Abstract. In this paper Gaussian models of retarded and accelerated anomalous 
diffusion are considered. Stochastic differential equations of fractional order driven by 
single or multiple fractional Gaussian noise terms are introduced to describe retarding 
. and accelerating subdiffusion and superdiffusion. Short and long time asymptotic 

£3 \ limits of the mean squared displacement of the stochastic processes associated with the 

solutions of these equations are studied. Specific cases of these equations are shown to 
provide possible descriptions of retarding or accelerating anomalous diffusion. 
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1. Introduction 

Anomalous diffusion occurs in many physical, chemical and biological systems [1-3]. 
In normal diffusion the mean squared displacement of diffusion particles varies linearly 
with time (Ax 2 (t)) ~ t. In some complex disordered media the diffusion becomes 
anomalous with (Ax 2 (t)} ~ t a , where the scaling exponent a ^ 1 characterizes the 
anomalous diffusion. For a < 1 the process is known as subdiffusion, and when a > 1 
it is called superdiffusion. Differential equations of fractional order are well-suited for 
describing fractal phenomena such as anomalous diffusion in complex disordered media. 
The constant memory and self-similar character of these phenomena can be taken into 
account by using the kinetic equations with fixed fractional order. 

There exist in certain locally heterogeneous media where diffusing processes do not 
satisfy the constant power-law type scaling behavior like anomalous diffusion. Such 
processes include the retarding and accelerating anomalous diffusion. Retardation of 
diffusion occurs in single-file diffusion where particles are constrained to move in single- 
file due to confined one-dimensional geometry, such as diffusion in zeolites [4], or in 
the anomalous diffusion that occurs on the biological cell membrane [5,6]. Possible 
causes for the presence of anomalous subdiffusion in biological systems include the 
presence of immobilized obstacles which hamper molecular motion by an excluded 
volume interaction, and the cytoplasmic crowding in living cells. On the other hand, 
membrane-bound proteins exhibit transition from subdiffusion at short time to normal or 
superdiffusion at long times [7] ; the diffusion of telomeres in the nucleus of mammalian 
cells shows accelerating subdiffusion [8]. Other examples are retarded and enhanced 
dopant diffusion in semiconductors [9-11], the accelerating superdiffusion of energetic 
charged particles across magnetic field in astrophysical plasma physics [12,13], and 
accelerated diffusion in Josephson junction [14]. Such processes have memory and fractal 
dimension that may vary with position, temperature, density, or internal parameters 
such as elasticity and viscosity. 

Anomalous diffusion having scaling exponent which varies with position and time 
was studied by Glimm and co-workers in the early 1990s [15, 16] in the multifractal 
modeling of heterogeneous geological systems. More systematic studies of transport 
phenomena with variable scaling exponents were carried out in the late 1990s and 
early 2000s. Stochastic processes with variable fractional order such as multifractional 
Brownian motion was introduced to model phenomena with variable memory or variable 
fractal dimension [17-19]. In order to describe systems with variable scaling exponents 
one may have to consider fractional differential equations of variable order. However, 
such variable order equations are in general mathematically intractable and can not be 
solved without numerical approximations [20-25]. 

There exists certain class of diffusion processes with non-unique scaling exponent 
that can be described by fractional differential equations of distributed order. The 
notion of distributed order differential operators was first introduced by Caputo in 
1969 [26]. A distributed-order fractional diffusion equation has its fractional order 
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derivatives integrated over the order of differentiation within a given range. Applications 
of distributed fractional order equations to fractional diffusion and fractional relaxation 
have been carried out by various authors [24-32]. 

This paper considers Gaussian models of retarded and accelerated anomalous 
diffusion. We introduce a class of multi-term fractional Langevin-type equations driven 
by single or multiple fractional Gaussian noise terms. These equations can be regarded 
as special cases of fractional Langevin equations of distributed order, and they can be 
used to describe retarded and accelerated anomalous diffusion. Detailed study of the 
short and long time asymptotic properties of the solutions to these equations are carried 
out. 



2. Multi-fractional stochastic differential equations 

In this section we introduce a class of multi-fractional Langevin-like equations of the 
following form: 

m n 

aiD ai x(i) = c i^v t eR,0 <a { <2 (1) 

i=l j=l 

where Cj > 0, D"' is the Riemann-Liouville or Caputo fractional derivative [33-37] which 
is defined for m — 1 < a < m as 

1 d m /■* f(u)du T . ... 

— / r — , Kiemann-Liouville 

Dff(t)={ r(m- a) d( T i„ (t -„)—« (2) 

' " ti)"*-"- 1 — /(ti)Ai, Caputo 



4l / <' " 

- «) Jo 



r(m -a) J K dt r ' 
The Gaussian noise (t) is defined by 

($»(*)) = o, (3) 

and 

= - s r 7j ' < 7i,7j < 2, (4) 



with 

^ = 2 8in(7r 7i /2)r(l- 7j -)' (5) 
If we let 7j • = 2 — 2^, where < iJj < 1 is the Hurst index associated with fractional 
Brownian motion, and dj = (2 sm(7rHj)T(2Hj — (t) can then be regarded as the 
derivative (in the sense of generalized functions) of fractional Brownian motion indexed 
by Hj. Note that the covariance of fractional Gaussian noise has the same algebraic 
sign as (2Hj — 1). For 1/2 < Hj < 1, the process exhibits long-range dependence with 
persistent positive covariance. On the other hand, when < H 3 , < 1/2, dj is negative 
and the process has anti-persistent correlation structure. For Hj = 1/2, or 7j = 1, 
it corresponds to white noise. Here we remark that \im Hj ^i/ 2 dj\t\ 2Hj ~ 2 = S(t) in the 
sense of generalized functions [38,39]. We thus see that there is no need to include in 
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the covariance of fractional Gaussian noise an extra term to cater for the white 
noise when H — 1/2 as given by [40] 

^ W*)6f») = MjHjQHj — — s| 2 ^- 2 + 4^|t - s\ 2H - l 5{t - s). (6) 

Here we would like to briefly discuss the covariance of fractional Gaussian noise 
[40-42] in terms of generalized functions. Just like for a proper description of fractional 
Gaussian noise, should not be defined pointwise for each t. Instead it needs to be 
considered the process as a linear functional in some test function space such as Schwarz 
space <^(R) of real- valued infinitely differentiable functions which decrease rapidly [38]. 
The generalized process / e ^(M), is a linear functional 

/oo 
ti(t)f(t)dt. (7) 
-oo 

£(/) is a generalized stationary Gaussian process with covariance given by the bilinear 
functional 

C(f,g) = U(m(9))= [°° [°° f(t)g(s)c(t-s)dtds 

J — oo J — oo 
/•oo 

/(t)dt / g(s)c(\t-s\)ds 



oo </ — oo 

oo POO 



/oo poo r _. 

f(t)dt c(s) g(t + s)+g(t-s) ds, (8) 
-oo J0 1 J 



where c(t — s) is a generalized function or distribution. For example, for H = 1/2, £ is 
white noise with c(t — s) = 5(£ — s). H ^ 1/2 corresponds to fractional Gaussian noise, 
with c(t), t > given by 

t 2H-2 f T^T^), 1/2 < If < 1 

1 J sin(^)r(2iJ - 1) 1 ni-^ (f) < // < 1/2 ' ^ ^ 

I, sin(7r/f) 

When 1/2 < if < 1 the covariance kernel in (8) can be regarded as the fractional 
integral of delta function, J 2H ^ 1 5(t) up to a multiplicative constant (2 sin(7r if) -1 . In 
this case, c(t — s) is locally integrable. It becomes fractional derivative of delta function 
for (2sin(7rii)- 1 D 1 - 2// 5(t) for < H < 1/2 since I a f = D~ a f [42]. 

Before we proceed further, a brief comment on the stochastic differential equation 
driven by fractional Gaussian noise will be given. Recall that stochastic calculus of Ito 
can not be used to define the integrals with respect to a stochastic process which is 
not a semimartingale. Fractional Brownian motion is not a semimartingale when the 
Hurst index of the process H ^ 1/2, or when it is not a Brownian motion. Due to its 
widespread application, the question on how to obtain a well-defined stochastic integral 
with respect to fractional Brownian motion has become a long standing problem which 
has attracted considerable attention [43,44]. Several methods which include Sokorohod- 
Stratonovich stochastic integrals, Malliavin calculus, and pathwise stochastic calculus 
have been proposed to overcome this difficulty (see references [44,45] for details). 
However, for application purposes, theory based on abstract integrals may encounter 
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difficulty in physical interpretations. As we shall restrict our discussion related to 
applications fractional Gaussian noise involving only persistent case 1/2 < Hj < 1 
(or < 7j < 1), the integrals with respect to fractional Brownian motion can thus be 
treated as the pathwise Riemann-Stieltjes integrals (see for example [45] and references 
given there). This allows one to handle such integrals in a similar way as ordinary 
integrals. 

Here we remark that the multi-term fractional order Langevin-like equation (1) 
can also be regarded as a special class of the following distributed-order fractional time 
stochastic equation: 



D v x(t) = ^(t), t > 0. 



(10) 



(11) 



with the distributed fractional derivative 

D v x(t) = [ (fi(a)Dfx(t)da, 
Jo 

where D" is the fractional derivative as defined by (2), the weight function (p(a), 
< a < 2, which satisfies ip(a) > is given by ip(a) = Y^T=i a $( a ~ Note 
that in general the weight function is a positive generalized function. The Gaussian 
noise £^(t) is the distributed-order fractional Gaussian noise defined by 

Z*(t) = [ UW(ri)di, (12) 

Jo 

with the weight function ^(7), < 7 < 2, which satisfies ^(7) > and is given by 
V'(t) = J2 f j=i c jfi(,7 ~~ Here we remark that in most of the examples considered 
subsequently, we shall restrict < 7 < 1. 

Solution of (1) can be solved formally by using Laplace transform method which 
gives 

A(s)x(t)-B(s)=i(s), (13) 
where x(s) is the Laplace transform of x(t) and 



A(s) = ^a iS a \ 



(14) 



and 



B(s) 



m [oti] 

^2 a i ^2 sk \ D RL kllx< ^) ' (Riemann-Liouville) 



k t =0 
foil 



i=l ki=0 



m 



(15) 



(0), 



(Caputo) 



Here \af\ denotes the largest integer smaller or equal to For Riemann-Liouville case, 



RL 



is the Riemann-Liouville derivative of order en,- — h — 1 evaluate at 



t=o 



t = 0. For the Caputo ( ki >(t) denotes fcjth derivative of x(t). For simplicity we 



assume the initial conditions x^ ki \0) = for all % 



, to, and 



D 



cti—ki — l 
RL 



X(t) 



t=0 



0, 
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such that B(s) = for both these cases. The Laplace transform of the Green function 
is then given by G(s) = l/A(s). Therefore, 

x(s) = G(s)i(s) = M. (16) 
The solution is then given by the inverse Laplace transform: 

x(t) = [ G(t - u)i(u)du. (17) 
Jo 

The covariance and variance of the process are given respectively by 

K{s,t)= f du f dvG(t-u)C(u-v)G(s-v), (18) 
Jo Jo 



and 



a 2 (t) = f f G(u)C(u - v)G(v)dudv 
Jo Jo 

= 2 G(u) / C(u — v)G(v)dudv. (19) 
Jo Jo 

Assuming t > s, (18) becomes 

K(s,t)= f dv\G(s-v)Gc(t-v) + G(t-v)G c (s-v)], (20) 
Jo L - 1 



with 



G c (t) = (G*C)(t) = [ G(t-u)C(u). 

Jo 



(21) 



Here we would like to remark that throughout this paper we assume the stochastic 
processes under consideration have zero means, so the mean squared displacement is 
equal to variance. These two terms will be used interchangeably in our subsequently 
discussion. 

Note that for the simplest case of (1) with m — 1, n — 1 and ol\ = a, 71 = 1, then 
for Riemann-Liouville (or Caputo) fractional derivative, 

D a x(t)= V (t), (22) 

where = rj(t) is white noise. For D^xif) = (or x(0) = 0), (22) defines 

t=o 

Riemann-Liouville fractional Brownian motion or type II fractional Brownian motion 
with Hurst index H, a — H + 1/2 [19]. The solution of (22) with the above boundary 
condition is given by 

x(t) = I a V (t) = [\t - uf-^du, (23) 

r(«) Jo 

with variance 

j.2a-l 4.2H 

a 2 (t) = j = s . (24) 

(2a-l)(r(a)) 2 2H{T{H+ 1/2)) 2 

Note that fractional Brownian motion of Riemann-Liouville type has a variance with 
the same time dependence as the standard fractional Brownian motion. In contrast to 
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the later, though it is self-similar but its increment process is not stationary. However 
it has the advantage that the process begins at time t — 0, and the Hurst index can 
take any value H > [19]. 

Another simple case is when m — 1, n — 1 and 71 = 7 7^ 1, which leads to a 
Gaussian non-stationary mono-fractal process with variance 

°" ^ = sin(77r/2)(2a - 7 )r(a)r(a - 7 + 1) ' ^ 
This process can be subdiffusion or superdiffusion, depends on whether 2a — 7 < 1 or 
2a - 7 > 1. 

In the subsequent sections we shall consider various specific cases of (1) for modeling 
accelerating and retarding anomalous diffusion. 



3. Accelerating Anomalous Diffusion 



One of the simplest model for accelerating diffusion can be obtained by using a special 
case of (1): 



""•'•in 

3=1 

The Green function is given by G(t) = and 



G 



c 



(*) = f 
Jo 



du 



(t-u 



,a-l 



r(a) 



^ ^2sin(7r 7 ,/2)r(l- 7j ) 



-j 

u <*-lj 



^ '2sm(7rV2)r(a-7i + l)' 
The covariance is given by 

K( S ,t)=^^( S ,t), 
J'=l 

and from (20) one gets 



Kj(s,t) 



2sin(7T7j/2) 



■1i 



(s-u)^ 1 (t-u) a - 
T(a) T(a - 7i + 1) 



■7j 



(t-w)"- 1 (s-u) Q - 
F(aj r(a - tJ + 1) 



2sin(7r 7i /2) 



r(a)r(a - 7j + 1) 7o 



g a— 7j+l^Qi— 1 



r(a)r(« - 7j + 1) 



Jo 



du(l-u) a ~ l (l - |w a " 7i ) 



"7j 



2sin(7r7,-/2) 



r(a + i)r(a- 7i + i) F(7j - " 1 + 



-7j+l^a— 1 



]»!> - 7 j + 2) 



F(l-a,l,2 + a-7j,s/t) 



(26) 



(27) 



(28) 



(29a) 



(29b) 



8 




<B> V 



Figure 1. The domain of different anomalous diffusion types 



where we have used [46], #9.111, page 1005. By using the following identity: 

r( 7 )r( 7 -a-/3) 



F(a,0,j,l) 



r(7-«)r( 7 -/3) : 



one gets the variance as 

n n 



n 

E 

j 



c 2 t 2a-v 



1 (2a - jj) sin(7T7 i /2)r(a)r(a - 7j 



(30) 



, (31) 



Note that the variance can also be obtained directly from (29a). 

For the case with n = 2, and 7i > 72, then the short-time and long-time limits of 
the mean-square displacement (or variance) are given by 

C 2 t 2«- 71 

At) 



sin(7r 7l /2)(2a - 7l )r(a)r(a - 7l + 1) 

c 2 t 2a- 72 



as t ->■ 



(32) 



as t — >• 00 



sin(7r 72 /2)(2a - 72 )r(a)r(a - 72 + 1) 

Therefore the process is accelerating subdiffusion if for i — 1, 2, < (2a — 7 «) < 1, and 
it becomes accelerating super diffusion when 1 < (2a — 7 j) < 2. For < (2a — 7 i) < 1 
and 1 < (2a — 72 ) < 3, the process begins as a subdiffusion and it accelerates to become 
a super diffusion. Figure 1 shows the various possible type of anomalous diffusion for 
different values of a and 7 . 

In order to have some idea about the type of process x(t) represents, let us consider 
some specific examples correspond to certain given values for a and 7 j. For (26) with 
a = 1, the ith component of covariance of the process associated with the solution 
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becomes 



Ki(s,t) 



2sin(7r 7i /2)r(2- 7i 



2-7* 



(t - u) 1 -^ + (t- uf-^ 
+ s 2 '^ - (t - sf-^ 



du 



2sin(7r 7i /2)r(3- 7i )r J (33) 

which is the covariance of the fractional Brownian motion (up to a multiplicative 
constant). Thus the covariance of the process K(s,t) = X^=i Ki{s, t) is the covariance 
of a mixed fractional Brownian motion (also called fractional mixed fractional Brownian 
motion by some authors) [47-49] which is the sum of n independent fractional Brownian 
motion 



x(t) = Y,*B Hi (t), 



(34) 



where — (2 — 7»)/2 is the Hurst index of the fractional Brownian motion B H .(t). 
Thus, in the mixed fractional Brownian motion model, anomalous diffusion begins with 
a lower diffusion rate can be represented by the fractional Brownian motion of lower 
Hurst index, and it is subsequently accelerated and is described by fractional Brownian 
motion of higher Hurst index. It is interesting to note if all are not equal to | (that is 
when they are all not Brownian motion), such a process is long-range dependent [50, 51]. 
The process satisfies a generalization of self-similar property called mixed self-similarity 
in the following sense: 



Y^B Hi {rt)±^Tr^B Hi {t). 



i=i 



The variance of the process is given by 



Just like the previous case, suffice to consider n 
time limits of the MSD are 



(35) 



(36) 



2. For if 2 > Hi, the short and long 



" 2 t 2Hl 



a 2 (t) 



sin(7rifi)r(2#i + 1) 



as t — > 0, 



as t — >• oo. 



(37) 



t sin(7r# 2 )r(2# 2 + 1) 
Thus the process behaves as accelerating superdiffusion (or subdiffusion) for 1/2 < Hj < 
1 (or < Hj < 1/2) with j = 1, 2. 

Here we would like to remark that there exists another process called step fractional 
Brownian motion [52-54] which can also be used to model both accelerating and 
retarding anomalous diffusion. Although the mixed fractional Brownian motion can 
only be used to describe accelerating anomalous diffusion, its mathematical structure is 
comparatively simpler than that of step fractional Brownian motion. 
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Next, we consider another special case with the exponents a and 7, satisfy the 
condition a — 7, = for a particular ith fractional Gaussian noise. (29a) now becomes 



c? 



KAs, t) = . ^ r— . / (s - u) Q_1 + (t - uf' 1 du 

2sin(7ra/2)r(a) J 1 



-,2 



s a + t a - (t - s) a ] . (38) 



2sin(7ra/2)r(a + 1) 

Since 7$ = 2 — 2if, = a, the process with the above covariance is fractional Brownian 
motion indexed by 2 — 2iJj. Note that the other components of covariance Kj(s,t), 
j ^ i are given by (29b) hence they are not fractional Brownian motion. The variance 
is given by 

n 



2sin(7ra/2)r(a + 1) 

+ V ^- f39) 

.^2shi(7r7 i /2)(2a-7i)r(a)r(a-7i + l)" 1 ^ 

In this case the type of anomalous diffusion depends on values of a and 2a— jj. Assuming 

< 7j < 1, then for < a < 1 and < 2a — 7,- < 1 one gets accelerating subdiffusion; 
when 1 < a < 2 and 1 < 2a — 7j < 2, the process is accelerating superdiffusion. The 
process accelerates (retards) from subdiffusion (superdiffusion) to become superdiffusion 
(subdiffusion) when 1 < a < 2 and < 2a — jj < 1 (or when < a < 1 and 

1 < 2a -7j- < 2). 

Finally, we note that if & is white noise, that is when 7^ = 1, then Ki(s,t) takes 
the following form: 

1 (a + 1)1 (a) 

which is just the covariance of "type II" or Riemann-Liouville fractional Brownian 
motion [19]. The variance for the process is given by 

o„2 f 2Q-l 

= (9 1VTV ^2 - ( 41 ) 
(2a- l)(r(a)) 2 

For n = 2, 71 = 1, and 72 < 1, one gets a process which is an accelerating subdiffusion 
ifl/2<o;<l;an accelerating superdiffusion if 1 < a < 3/2, and finally it represents 
a process which accelerates from a subdiffusion to a superdiffusion if l/2ct < 1 and 
72 < 2a - I. 

From the above discussion it is noted that a simple fractional Langevin equation 
driven by a single fractional Gaussian noise results in an anomalous diffusion. However, 
when the process is driven by more than single fractional Gaussian noise term, the 
resulting process can be accelerating subdiffusion or superdiffusion, depending on the 
fractional order of the noise terms and the derivative term. Thus, the interplay between 
multiple driving fractional Gaussian noise terms in the fractional Langevin-like equation 
leads to a simple Gaussian model of accelerating anomalous diffusion. 
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4. Retarding Anomalous Diffusion 



There exists another type of anomalous diffusion that slows down with time. In other 
words, for short times the mean square displacement of such a diffusion process varies 
as t ai , and it varies as t° 2 with a 2 < ol\ for long times. Retarding anomalous diffusion 
occurs in various physical and biological systems. For examples, in the single-file 
diffusion that occurs in cell membranes and narrow channels [55,56], and the retarding 
anomalous diffusion in semicondutors [11,57]. 

A simple Gaussian model for retarding anomalous diffusion can be obtained by 
considering a special case of (1) which takes the form of following fractional stochastic 
differential equation: 



ai 



D ai x{t) + a 2 D a2 x(t) = CjZ yj (t), < «2 < Oil < 2. 



(42) 



h=i 



For both Riemann-Liouville and Caputo case, the Laplace transform of (42) gives 



diS 



ai 



a 2 s 



x(s)-B(s) = J2^M> ( 43 ) 

3=1 

where B(s) is given by (15). For simplicity, we again choose the initial conditions such 
that B(s) = 0. The Green function is then given by the inverse Laplace transform of 



G(s) 



ai-s ai + a 2 s° 2 ai s Q i- a 2 + (a 2 /ai) 
such that its inverse Laplace transform gives 



G{t) t E ai — a2tai 



(44) 



(45) 



where 



3=0 



li > 0, v > 0, 



(46) 



is the Mittag-Leffler function 

En o t 'i j. 
j= iC jW z^ one gets 



Using (20) and (21), and substituting C(t) 



^ io «i V °i / 



3 

n c 2 

ai 



ai— 02,0:1 — 7j+l 



di 



The covariance K(s,t) = YTj=i c j^j( s ^) ^ s gi ven by (20) with 



Kj(s,t) 



du 



-E, 

ai 

(t - u) ai -^ 
ai 



01 — 02,01 



° 2 ( V- 

(s-u) 

di 



E ai — ct2,ai— -fj+1 



(t-u) 

di 



r(i - 7i ) 

(47) 
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+ 



du 



(t-u) ai - 



-E, 



(s - u) ai ~^ 



ct\—a.2,oti 



0-2 / 



U 



,«l-«2 



E, 



ai— Q2,Qi— 



(48) 



The variance is given by cr 2 (t) = YTj=i C j-^.?'(M) 

'(t-u) 01 -^- 1 



K 



(*,*)= 2 r 

Jo 



du 



E, 



E 



ai— 02,01— 7j+l 



Q>2 / 



u )ai-aa 

- w) ™ 



(49) 



Note that the covariance and variance given above can not be evaluated. However, 
by using the following asymptotic properties of Mittag-Leffler function [58], 



N 



and 



-E 



T(u — nji) 



+ 0[\z 



-1-N 



arg(z)|< (i-|)tt, 



z — > oo, 



r» 



+ 0(z), ^^o, 



(50) 



(51) 



it is possible to obtain the short and long time behavior of the variance. In the case of 
single fractional Gaussian noise, one has for 7 < 2a 2 , 

f 2ai- 7 t ^Q 



a\t) 



J.202— 7 



t — y 00 



(52) 



Since a\ > a 2 , the process is a retarding subdiffusion (or superdiffusion) if for % = 1,2, 
< a{ — 7/2 < 1/2 (or 1/2 < ccj — 7/2 < 1). In the case when there are more than one 
noise, the dominant terms for the short time limit and long time limit for the variance 
are ~ ^mm(2ai-7j) anc i ^ f max(2o; 2 —yj ) reS p ec ti V ely Thus we see that the double-order 
fractional stochastic equation driven by single fractional Gaussian noise can be used to 
model retarding anomalous diffusion. The lower order fractional derivative term in (42) 
plays the role of a damping term which slows down the diffusion. 

The main disadvantage of using (42) for modeling retarding anomalous diffusion 
is that in general the covariance and variance of the underlying process can not be 
calculated explicitly. We would like to find a particular case of (42) such that its solution 
is a process with covariance and variance that can be completely determined. For this 
purpose we consider the following double-order fractional Langevin-like equation: 

ctiD ai x(t) + a 2 D a2 x(t) = Cl ^(t) + c 2 ^ 2 (t), < a 2 < an < 2. (53) 

The two independent fractional Gaussian noises and C, 2 (t) are chosen such that for 
t > s they have zero mean and the following covariance 

(t - s)"-* -1 



-Sij = Ci(t - s)8ij, ij = 1,2. 



(54) 



We remark that the fractional Gaussian noise with covariance given by (54) is selected 
based on practical purposes as it gives the required results as well as provides a more 
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manageable solution to (53). As a result (54), one can verify that the Laplace transforms 
of the covariance of the fractional Gaussian noise C(t) = c\Ci(t)+c\C 2 (t) and the Green 
function G(t) satisfy the following relation: 

G(s)C(s) = s- u . (55) 

The inverse Laplace transform of (55) is 



G(t) * C(t) 



(56) 



The covariance of the process associated with the solution of (53) can be calculated 
using (20) and (56): 



K(s,t) 



du 



o 



G(t-u) - \ ^G(s-u) K 



r(iz) 



— du 
Jo 



(t-u) 



ai — 1 



-E, 



.^(t_ u )«i-« a 



u 



,J/-1 



r(iz) 



(s-u) 



-E, 



ai 



a 1 -a 2 ,a 1 



a 2 , 
— s 

Its variance is given by the following series expansion: 



(t-u) 



a\t) = 2 f 
Jo 



du 



aiT(is) 



(t-u) 



ai+U—2 rn 



( ~°2 A +n 



-^(t-u) ai - a2 



(57) 



(Oi — Q2)+«l+l y — 1 



a iT(^) ^ (n(ai — a 2 ) + ol\ + v — l)r(n(a;i - a 2 ) + «i) 
From (58) one gets the short and long time limits of the variance as 

-t ai+v -\ t^O 



(58) 



a\t) 



a 1 T(iy)T(a 1 )T(a 1 + iy-l) 



t a *+"-\ t^oo 



(59) 



a 1 T(u)T(a 2 )T(a 2 + v - 1) 

The process with covariance (58) and variance (59) is a retarding subdiffusion or 
superdiffusion depending on the on the values of ctj + v — 1, % — 1, 2. Thus, one can 
regard the term a 2 D a2 x(t) in (53) as a damping term which slows down the anomalous 
diffusion. 

A special case for which the covariance has a closed form is when v — 1. For t > s, 
the covariance (57) becomes 



K(t,s) 



du 



—ct-ur^E, 

ai 
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(60) 
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It is interesting to note that the covariance consists of three terms of Mittag-Lemer 
functions of same order, and the time variables t and s enter the covariance expression 
(60) in a form similar to that of fractional Brownian motion. Thus it is not a coincidence 
that the asymptotic short and long time limits of the covariance are given by 

1 



K(t,s) 



aiT(ai 



f a 2 + g a 2 



t,s ->■ 

t, s, \t — s\ — ?■ oo 



• (61) 



a 2 T(a 2 + 1) 

Note that these are just the covariance of fractional Brownian motion indexed 
respectively by a±/2 = 1 — Hi and 0:2/2 = 1 — H 2 , with < ifj < 1, % — 1, 2. Since both 
the short and long time limits are fractional Brownian motion, one can conclude that 
the stochastic process in this case is long-range dependent except when cci = a 2 = 1/2 
or when both the limiting processes are Brownian motion [50, 51]. 
The variance is given by 



<J 2 {t) = -t ai E ai . 



ai 



cxi —0:2,0:1+1 



a 2 
— i 
ai 



(62) 



The long and short time limits are given by 



CLiYioii + 1) 

2 



t c 



-t c 



t ->■ 
t — > 00 



(63) 



a 2 T(a 2 + 1) 

From the above results one has for < a 2 < ot\ < 1 a subdiffusion process which slows 
down with time, or a retarding subdiffusion. On the other hand, if 1 < a 2 < ot\ < 2, 
the process is a retarding superdiffusion. 

It would be interesting to see whether (53) can be used to describe accelerating 
anomalous diffusion as well. Suppose we replace the condition (55) by the following: 

G(s)C(s) = s~ v + s~ K , (64) 

such that 

C(s) = ais ai ~ u + ais ai - K a 2 s a2 - u + a 2 s a2 ~ K . (65) 
Inverse Laplace transform of (65) gives 



C(t) = ai 



^v—ai — l ^-k— 01 — 1 

— + 



Otl) 



+ a 2 



-j-l/—a.2 — l ^-K — — 1 

+ 



T(u - a 2 ) T(k - at 2 ) 



(66) 



Y(v — ot\) T(k 
The variance of the resulting process is 

ah) = al(t) + a 2 K (t), (67) 

with <rl(t) given by (58) and similarly for cr%(t) with v replaced by k. Thus for v < k, 
one gets 

2 



Oir(i/)r(o;i)r(ai + v — 1) 
2 

a 2 r(K)r(o;2)r(a;2 + k—1) 



(68) 
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It is interesting to note that the process with the above variance represents accelerating 
subdiffusion if a-y + v < a 2 + k, or k > a.\ — a 2 + v > 2a\ — a 2 . However, to achieve such 
an accelerating subdiffusion it is necessary to consider four fractional Gaussian noise 
terms in (53), a situation that may be difficult to realize in practice. 

On the other hand, we recall that for (42) with more than one fractional Gaussian 
noise, the dominant term of the variance for the associated process in the short and long 
time limit is respectively varies as t min ( 2a i-ij) an d t max ( 2a 2--yi) _ jf we consider the case 
with cki > a 2 and 71 > 72, one then has 2a± — 72 > 2cci — 71 and 2a 2 — 72 > 2a 2 — 71. 
As a result, 



Thus, it is possible to use (53) to model accelerating anomalous diffusion provided 
2a 2 — 72 > 2ai — 71 or a.\ — a 2 < (7! — 72V2. In other words, the fractional stochastic 
equation (53) can be used to model retarding and accelerating anomalous diffusion by 
appropriate choice of the order of the fractional derivatives and fractional Gaussian noise 
terms. Note that retarding diffusion such as single-file diffusion can also be modeled by 
fractional generalized langevin equation [59]. 

5. Concluding remarks 

We have shown that it is possible to model both accelerating and retarding anomalous 
diffusion by using fractional Langevin-like stochastic differential equations driven by one 
or more terms of fractional Gaussian noise. The solutions associated with some specific 
cases of these equations turn out to be some interesting processes in the short and long 
time limits. For example, two types of fractional Brownian motion, namely the usual 
standard fractional Brownian motion and the Riemann-Liouville fractional Brownian 
motion are the asymptotic processes of special cases of the model. This model also 
includes another interesting process, namely the mixed fractional Brownian motion, 
which is a simple process for describing accelerating sub- and super-diffusion. 

We note that the stochastic differential equations in our model can be regarded as 
fractional Langevin-like equation of distributed order (10) with weight function consists 
of delta functions. One may want to consider cases with different type of weight functions 
in (10), such as uniform or power-law weight functions. However, from the results of 
our previous study on fractional Langevin equations of distributed order with uniform 
and power-law type of weight functions indicates that such equations in general do not 
have closed solutions even for the case of simple fractional Langevin of distributed order 
driven by white noise [32]. One thus expects the situation to be even more complex when 
weight functions other than the delta functions are used for the multi-term fractional 
Langevin equation with more than one fractional noise terms. 

One question of interest is that whether it is possible to model accelerating and 
retarding anomalous diffusion based on (26) and (42), using different Gaussian noise. If 




(69) 
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one is only interested in the asymptotic limits of the mean squared displacement of the 
stochastic process, then instead of using fractional Gaussian noise in the stochastic 
differential equations (26) and (42), Gaussian noise with covariance which has the 
"correct" asymptotic limits of power-law type can be used. For example, Gaussian noise 
with covariance which varies as At^~ 2 / (1 + At^~ u / B), A and B are positive constants, 
< pi, v < 2. Such a covariance has respectively short and long time limit At^~ 2 and 
Bt v ~ 2 respectively Another example is the noise of Mittag-Lemer type with covariance 
that of the form At^~ 2 E^ v ^ + i (—At^ u /B) ) fi > v, which has short and long time 
asymptotic limit At^ 2 /T(fi + 1) and Bt u ~ 2 /T(v + 1) respectively. These two examples 
give the possible alternatives to the fractional Gaussian noise £ 7j as the driving noise. 
More discussion related to these cases will be given in a forthcoming paper [60]. 

The usual characterization of anomalous diffusion uses its mean squared 
displacement (or variance in the context of this paper). However, it is well-known that 
even for a Gaussian model mean squared displacement does not determine completely 
the underlying stochastic process and hence the mechanism of the anomalous diffusion. 
Recent advances in particle tracking devices allow experiments to track the trajectories 
of single molecule or nanoparticle in complex systems such as cells in a biological 
system. Mean squared displacement obtained from the time series data gives the scaling 
exponent of the anomalous diffusion undertaken by such particles. Comparison of the 
experimental data so obtained with various models of anomalous diffusion allows one to 
distinguish the different possible subdiffusion mechanisms. In particular, information 
of single-particle trajectories allows one to test the validity of ergodic property of the 
associated diffusion process. A stochastic process is said to be ergodic if the ensemble 
average of certain physical quantity such as mean squared displacement measured in bulk 
coincides with the time average of the same quantity over sufficiently long time from 
the single-molecule time series. Examples of ergodic process are Brownian motion and 
the standard fractional Brownian motion. Another process of interest which is ergodic 
is mixed fractional Brownian motion which is the sum of two independent fractional 
Brownian motion. On the other hand, Riemann-Liouville fractional Brownian motion 
and heavy tailed continuous-time random walk are non-ergodic [61-64]. 

By using the single-molecule data the comparisons of experimental data based 
on various models of anomalous diffusion such as the continuous-time random walk, 
fractional Brownian motion, fractional Levy stable motion, etc. have been carried out 
by various authors recently [65-72]. For examples, diffusion of beads in entangled 
F-actin networks, diffusion of at shorter times exhibits continuous-time random walk 
behaviour [65]. However, the analysis of the data of the anomalous diffusion in crowded 
intracellular fluid such as cytoplasm of living cells rules out continuous-time random 
walk and favours fractional Brownian motion [66]. Analysis of single particle tracking 
data of lipid granules in yeast cells by Tejedor et al [67] seems to rule out continuous-time 
random walk and shows agreement with fractional Brownian motion; but a subsequent 
study [68] shows that at short times the granules perform continuous-time random walk 
subdiffussion while at longer times the motion is consistent with fractional Brownian 
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motion. Various analyses of the biological data describing the motion of individual 
fluorescently labeled mRNA molecules inside live E. coli cells, a well-known experiment 
first conducted by Golding and Cox [69] do not lead to consensus result. Magdziarz 
et al [70] show that fractional Brownian motion as the underlying stochastic process; 
but subsequently Burnecki and Weron claim that the data follows fractional Levy 
stable motion [71]. According to Kepten et al the experiments on telomeres in the 
nucleus of the mammalian cell exhibit fractional Brownian motion [72]. Recent study 
by A.V.Weigel et al [73] on the physical mechanism underlying Kv2.1 voltage gated 
potassium channel anomalous dynamics using single-molecule tracking shows that both 
ergodic (diffusion on a fractal) and nonergodic (continuous-time random walk) processes 
coexist in the plasma membrane. Though it is widely recognized that the diffusion 
pattern of membrane protein displays anomalous subdiffusion, however, there is still no 
agreement on the mechanisms responsible for this transport behaviour. Currently there 
is still no consensus on whether heavy tailed continuous time random walk, fractional 
Brownian motion, fractional Levy stable motion or some other stochastic processes 
can provide the correct description to anomalous diffusion in some biological systems. 
Thus, it is important to make use of data and information other than the mean squared 
displacement or the anomalous diffusion exponent to determine the type of mechanism 
and the stochastic process describing the anomalous diffusion. We hope that some of 
the processes considered in this paper may be of relevance in describing the anomalous 
diffusion in biological systems. 

Finally we remark that it would be interesting to investigate the Fokker-Planck 
equations associated with the processes considered in this paper. The mean first passage 
time for some of the simpler cases can also be studied [60] . 
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